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ABSTRACT 



A numerical code, USPOTF2, has been formulated to solve for the potential flow 
for two airfoils executing unsteady motions in an inviscid incompressible flow medium. 
This code is an extension of an existing code U2DI1F, which does the same calculations 
for the single airfoil case. The technique uses the well known Panel Methods for steady 
flow and extends it to unsteady flow by introducing a wake model which creates a non- 
linear problem due to the continuous shedding of vortices into the trailing wake. The 
presence of the second airfoil introduces a set of non-linear coupled equations for the 
Kutta condition. Numerous case-runs are presented to illustrate the capability of the 
code. The case of the step change in angle of attack is compared with Giesing’s work. 
All other case-runs are illustrated together with the results for the single airfoil case. 
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THESIS DISCLAIMER 

The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every' effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. INTRODUCTION 



A. GENERAL 

In this paper , a numerical method is formulated to solve for the flow about two 
two-dimensional airfoils which are arbitrary located and are performing an arbitrary time 
dependent motion in an inviscid incompressible fluid. The original work by Teng [Ref. 
1] is for a single arbitrarily defined airfoil in the same potential flow condition. The ex- 
tension of Teng's code to two bodies is considered here. Where possible, the same con- 
vention and notation are adopted. As for the single airfoil case, all velocities are 
non-dimensionalised with respect to the free stream and all lengths with respect to the 
chord length. 

A new subroutine (SUBROUTINE NEW'P OS) is added to transform either of the 
two local coordinate systems to the global coordinate system. A modification is entered 
into the treatment of the Kutta condition to make it consistent with the unsteady flow 
Kutta condition treatment. Subroutine COFISH is deleted and its role is taken over by 
Subroutine COEF which normally performs the formulation for the flow tangency con- 
dition for the unsteady case. A more accurate method is also introduced to obtain the 
velocity potential. The reader is referred to the work of Kramer (Ref. 2] for further 
improvement of the original code. 

This documentation is set up as for the original documentation in order that it will 
be easier to follow and cross-referenced. W'hile it is the intent of the author to keep this 
thesis as complete as possible, the reader is well advised to review Reference 1 for the 
work involved in the single airfoil case; no special effort will be made to reproduce it 
here. 

B. BASIC THEORY AND APPROACH 
1. Steady Flow Problem 

The treatment for the two airfoils case in steady incompressible flow follows 
closely to the single airfoil case. The governing equation, as for the single airfoil case , 
follows from the Conservation of Mass and the Condition of Irrotational Flow. 

The continuity equation of an incompressible fluid (div V = 0) and the condi- 
tion of irrotational flow ( curl V = 0) leads to the well known Laplace Equation. 

div( grad <? ) = 0 (1-1) 
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with <t> denoting the disturbance potential for the velocity. This is seen to be the classic 
Neumann problem of potential theory with the usual problem of defining the boundary 
conditions. The boundary conditions for the disturbance potential <£ are that its gradi- 
ent normal to the surface be equal to the normal velocity of the surface and that its 
gradient vanish at infinity; that is 

V0 • / 1 =V r »rtonS (1.2) 

lim V</>( P ) -*• 0 (1.3) 

P-> oo 

where V r is the resultant velocity of a point on the body as seen from the inertia frame 
of reference, n is the outward unit vector normal to the body, S is the body surface and 
P represents a general point. Equation 1.2 holds for both airfoils i.e. on both surfaces. 

The pressure coefficient is obtained through the Bernoulli's equation which de- 
rives from the Momentum equation. 

The approach adopted is associated with Hess and Smith [Ref. 3] who devised 
the popularly known PAS EL method in the early sixties. In words, the boundary or 
airfoil surfaces S, and S 2 about which the flow is to be computed is approximated by a 
large number of surface elements whose characteristic dimensions are small compared 
to those of the body. Over each surface element, a uniform source distribution and a 
uniform vorticity distribution is placed. The source strength (q.) varies from element to 
element, while the vortex strength (y,) is the same for all elements in the same airfoil but 
is different across the airfoil. The singularity strengths are determined from the flow 
tangency condition on both body surfaces and the two Kutta conditions at both trailing 
edges. With the determination of the singularity strengths, the relevant aerodynamic 
data can then be subsequently computed. 

2. Unsteady Flow Problem 

The unsteady problem is similar to the steady flow problem in that they both 
have the same governing equation viz. the Laplace equation, and that for both problems, 
the pressure and velocity are decoupled so that the velocity and pressure calculation can 
be computed separately and consecutively. 

This problem differs from the steady flow in that another model is required to 
simulate the continuous shedding of vorticity into the trailing wake. The existence of a 
vortex sheet behind the airfoil can be explained by the Helmholtz theorem which is 
basically a statement of the Conservation of Vorticity. This requires that any change in 
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the circulation around the airfoil must be matched by an equal and opposite vortex 
somewhere in the fiowfield. The presence of the countervortices provides the flow with 
a kind of a memory in that the flow at a particular time is affected by the bound circu- 
lation of the past. It is this non-linearity that distinguishes the numerical technique from 
the simple steady flow problem of solving N linear equations in N unknowns. 

The solution technique requires an iterative type solution. The present ap- 
proach follows closely the original panel method of Hess and Smith as described in the 
steady flow development, while with regard to the modelling of the wake, it adopts the 
procedure advocated by Basu and Hancock [Ref. 4). A uniform source distribution ( q) k 
and a uniform vorticity distribution [>’(/)]« as for steady flow is placed on each panel at 
time r* where j denotes the panel number and / the airfoil number. The wake consists 
of a single vorticity panel attached as an additional element on each airfoil through 
which the vorticities are shed into the respective wake and a series of point vortices 
which are being convected downstream with the fluid. A uniform vorticity distribution 
of strength (■/„(/))* is placed on the wake panel of each airfoil. This panel is further 
characterised by its length A(/) 4 and its inclination ( with respect to the respective 

local frame of reference. After each time step, the vorticity of the wake panel is con- 
centrated into a single point vortex and convected downstream. Simultaneously a new 
wake panel is formed. The downstream wake of point vortices is thus formed by the 
shed vorticity of previous time steps. 

C. SCOPE 

Chapter 11 extends the original code to handle the steady flow problem for two 
airfoils set at different relative distances and angles of attack. 

Chapter III deals with the unsteady problem for the two airfoils system. It intro- 
duces a new subroutine and a more accurate method of calculating the velocity poten- 
tial. The Kutta condition for the two airfoils system is specially treated as its unique 
problem of a non-linear coupled system is not seen in the single airfoil case. 

Chapter IV describes the computer program, its essential capabilities and limita- 
tions, its associated subroutines, its input requirements and its associated output print- 
out. 

Chapter V presents the results of some case-runs. Of interest is a comparison case 
of a step change in angle of attack (AOA) with Giesing [Ref. 5] for the same airfoil 
undergoing an impulsive start at the same AOA. Case-runs will also be run with both 
airfoils at large distances apart to compare with the single airfoil case. In addition, ex- 
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ample cases for which no comparisons exist are given, to exhibit the capability of the 
method. 

Finally, Chapter VI concludes with future development efforts and the application 
potential of this numerical method. 
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II. STEADY FLOW PROBLEM FORMULATION FOR TWO AIRFOILS 

A. GENERAL 

The modification work for the steady flow is straightforward. The revised program 
allows for two arbitrarily defined airfoils placed at an arbitrary distance set at different 
angles of attack. For reason of simplicity, the number of panels and nodes and the pivot 
location are set to be the same for both airfoils. 

B. FRAMES OF REFERENCE 

Three frames of reference are involved in the two two-dimensional airfoils' case in 
steady flow. These are three inertia frames of references as indicated in Figure 1. 

The first inertia frame of reference (also known as global frame of reference) is set 
at the pivot position of the first airfoil with the X-axis pointing in the direction of the 
free-stream velocity. The two other inertia frames of references, henceforth, will be 
known as two frozen local frames of referencel (jc,y, and x 2 y 2 ) are fixed respectively to 
each airfoil with the x-a\is coinciding with the chord line originating from the respective 
leading edge. The two local frames of reference are set apart by XShift and YShift on 
the global frame of reference. In steady flow, the fluid velocity and pressure depend only 
on the spatial coordinates (X.Y) and not on time. 

The two airfoils are defined in the local coordinate system as input data for sim- 
plicity and are then transformed to the global coordinate system through a knowledge 
of the relative positions of the 2 airfoils' pivots positions and the respective local angles 
of attack. 

C. STEADY FLOW PANEL METHODS 
1. Definition of nodes and panels 

Each airfoil surface is divided into n straight line segments called panels by 
(n-i-1) arbitrarily chosen points called nodes. The numbering sequence begins with 
panel 1 on the lower surface at the first airfoil trailing edge and proceeds clockwise 
around the airfoil contour so that the last panel on airfoil 1 ends on the upper surface 
at the trailing edge. This numbering sequence then proceeds in a similar fashion for the 
second airfoil ( See Figure 2). As with the single airfoil case, the numbering sequence 

1 For the steady case, the notation 'frozen' will be dropped 
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